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The many body quantum dynamics of dipolar coupled nuclear spins / = 1/2 on an otherwise 
isolated cubic lattice are studied with nuclear magnetic resonance (NMR). By increasing the signal- 
to-noise ratio by two orders of magnitude compared with previous reports for the free induction 
decay (FID) of in CaF2 we obtain new insight into its long-time behavior. We confirm that the 
tail of the FID is an exponentially decaying cosine, but our measurements reveal a second universal 
decay mode with comparable frequency but twice the decay constant. This result is in agreement 
with a recent theoretical prediction for the FID in terms of eigenvalues for the time evolution of 
chaotic many-body quantum systems. 
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One of the simplest NMR experiments concerns the 
flipping and recording of the nuclear magnetization, i.e., 
the measurement of the free induction decay (FID). Quite 
to the contrary, rigorous theory to calculate this decay 
is often lacking since large numbers of nuclear spins are 
interacting, similar to the situation in electronic mag- 
netism with exchange coupled electronic spins. Often, 
one can only estimate the short-time behavior of the FID, 
but is unable to find an analytical expression for the en- 
tire decay. If the FID is limited by the life-time of the 
nuclear levels due to a coupling to a thermal bath, the 
" lattice" , the decay is simply exponential and given by 
the spin-lattice relaxation time Ti. This can be the case 
in liquids where rapid motion averages, e.g., the inter- 
nuclear magnetic dipole-dipole interaction, but, at the 
same time, couples the nuclear spins to the large thermal 
bath of motional degrees of freedom. In solids, where 
the spin-lattice relaxation time is often rather long, the 
magnetic dipole interaction leads to a rapid decay of the 
FID leaving the spin system far away from thermal equi- 
librium for times of the order of Ti. In this case, where 
the FID or other spin coherences decay according to the 
time-evolution of the acting Hamiltonian, one expects on 
general grounds the decay not to be exponential. How- 
ever, that is what is often observed with experiments 

An ideal system for the investigation of dipole cou- 
pled nuclear spins is CaF2 where spin 1 /2 fluorine nuclei 
are located on a simple cubic lattice (the low abundance, 
small moment Ca spins can be neglected) . Clean crystals 
are easily available with a nuclear spin-lattice relaxation 
time Ti of minutes, which practically leaves only one rel- 
evant time scale set by the dipolar coupling (~ 20iis). 
Given the simple nature of the material and the challeng- 
ing physics, there has been persistent interest in CaF2 
since the early days of magnetic resonance [iHSl [MS] . 

While most experimental and theoretical work fo- 
cussed on the short-time behavior of the dynamics 
[2 [51 110H12| , more recent research concerns the long-time 



behavior of NMR signals [T3lfT6] . In particular in re- 
cent work of Fine [15l [1^ , the long-time behavior of the 
system is addressed based on the notion of microscopic 
quantum chaos. 

The dynamics of a lattice of classical spins can be de- 
scribed by a set of angles {(j>i,Oi} where the i-th pair of 
coordinates describes the orientation of the i-th spin with 
its tip on the surface of a sphere. In case of a dipole cou- 
pled system the equations of motion controlling the time 
evolution of the system are nonlinear. For a large num- 
ber of spins this eventually leads to a phase space that 
is dominated by chaotic regions. While a Markovian de- 
scription of such a system is usually applicable only for 
times much larger than the mean free time, this is not 
true for ensemble averaged quantities |17j . An ensemble 
of spins can therefore be described in terms of Brownian 
motion as soon as the spin system has lost memory of its 
initial configuration. Thus, the long-time behavior of the 
ensemble averaged quantities can be obtained by solving 
a correlated diffusion equation on a spherical surface |15) . 

Recently, this concept has been generalized to quan- 
tum spins by Fine [151 116j , and he predicted a universal 
long-time behavior of FID and the decay after a spin 
echo. 

For a system of spins 1 /2 on any Bravais lattice with 
the dipolar Hamiltonian 
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in a strong magnetic field (so that only the secular part 
has to be retained) the long-time behavior is given by 
G {t) ~ e"'''"* cos (w^i + (/ict)- This reduces to a 

single mode after sufficiently long times, i.e.. 
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In this conjecture, the decay modes in G{t) follow from 
a set of solutions of a correlated surface diffusion equa- 
tion in the finite parameter space of single spin variables: 
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fa{t,x) = e~'*"'*MA<,(a;), where Xa = 7o- + iwo- is an in gen- 
eral complex eigenvalue of the integro-differential opera- 
tor associated with the modified diffusion equation and 
u\{x) its corresponding eigenfunction. In particular, the 
conjecture proposes that does not depend on the ini- 
tial configuration of the spin system. It was argued by 
Fine [TS] that due to the Markovian nature of ensemble 
averages [17] one expects not only the decay constant of 
the slowest mode (cr = 1) to be of the order of the inverse 
decay time, i.e., 71 ~ 1/r, but also the difference to the 
second slowest exponent should be of the order of l/r, 
i.e., 72 — 71 ~ 71. There are no predictions for 102- 

The early work of Engelsberg and Lowe [5] is in rea- 
sonable agreement with this theory, but inconclusive with 
regard to these predictions that meanwhile triggered new 
experiments, e.g., with hyperpolarized polycrystalline 
solid xenon for better signal-to- noise ratio [4]. These ex- 
periments, as well as recent ones for CaF2 show that 
the long-time behavior is indeed universal since it does 
not depend on various initial preparations of the spin sys- 
tem, in agreement with Fine's theory. However, only a 
single mode Ai in the exponential decay was found. 

Here, we report on new experiments on CaF2 that, due 
to the increase in signal-to-noise ratio by two orders of 
magnitude over previously reported experiments, reveal a 
second decay mode in agreement with Fine's predictions 
[inillSI- Therefore, our results favor his theory over other 
theories using a memory function approach T^, since 
these predict only a single mode. 

The 5 X 5 X 10 mm'^ CaF2 crystal was obtained from 
Mateck, Germany. Impurities are stated to be below 
2 ppm. The crystal was oriented by X-ray diffraction. 
A home-built NMR probe was used in order to align the 
crystal's axes (110) or (100) parallel to the applied mag- 
netic field (7.06 T). The probe was set to a resonance 
frequency of 283.383 MHz and operated at 20 K. The 
quality factor Q of the resonance circuit was 240, the 7r/2 
pulse length was 5 jis. The Ti at 20 K was determined to 
be 76 s. 

With typically 80 scans our maximum signal to noise 
was about 10^. In order to record the signal, various 
attenuations were used to prevent the preamp from sat- 
urating and to allow for an appropriate load of the 16 bit 
digitizer. Additionally, the first 20iis where measured 
using a spin locking sequence [TB]. Finally, the tail of 
the FID was measured with the variable attenuator set 
to zero but still infront of the preamp to minimize addi- 
tional phase shifts. 

In order to determine and remove the influence of the 
inhomogeneous magnetic field on the decay, we measured 
a ^^F FID on a CqFq sample (intrinsic T2 ^ 1 ms) of 
about the same volume and geometry as the crystal. On 
the time-scale of the CaF2 FID the ^^F decay could be 
approximated by a Gaussian with a standard deviation 
of 632(3) jis. The CaF2 data were hence multiplied with 
a rising Gaussian function to remove the inhomogeneous 
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FIG. 1. (a) Free Induction Decay of ^^F in a CaF2 single crys- 
tal (blue line) for B parallel (110) at 7.06 T. Note that the real 
part of the FID changes sign at each of the minima in the log- 
arithmic plot (b). The first 20iis of the FID where measured 
using a spin-locking technique, c.f. [18]. The time intervals 20 
to 57 ps and 58 to 87 ps where measured at 20 K with 25 and 
12 dB attenuation in front of the preamp, respectively. The 
tail of the FID was measured without attenuation. The sig- 
nal was corrected for the inhomogeneous decay as measured 
by CeFe as well as for a small residual offset (700 Hz) that 
was precisely determined from the time dependence of the 
signal's phase. The data is in good agreement with Engels- 
berg and Lowe's measurements but has nearly three orders 
of magnitude better signal-to-noise ratio. Due to the higher 
signal-to-noise ratio, we were able to fit the FID to Eq. ^ in 
the time frame 160 to 250 11s (shaded in gray) . 



broadening. Since the linewidth hampers a precise de- 
termination of the exact resonance frequency by means 
of a simple Fourier transform, the actual offset was de- 
termined by analyzing the time-dependence of the phase 
(/)(t) of the complex valued FID. The offset of 700 Hz was 
removed from the FID. Phase changes due to the variable 
attenuator were taken into account. 

The resulting signal is shown in Fig. [l] linear (a) and 
logarithmic scale (b) . The time origin of the FID was as- 
signed to the center of the 7r/2 pulse. The large dynamic 
range of the FID is only seen in the logarithmic plot in 
Fig. [ijb), where we show the magnitude of the NMR 
signal's real part. 

We find good agreement with Engelsberg and Lowe 
for times accessible within their experiment, i.e., up to 
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FIG. 2. Extraction of two decay modes from the FID for B parallel to (110). (a) After subtracting the fit or first Mode /i as 
determined by fitting the tail of the FID (c.f. Fig. [T]) to Eq. ([2|, a second mode /2 becomes apparent. This mode is replotted 
in (b) along with a fit according to Eq. ([2| in the time intervall 76 to 136 jis (shaded in gray). The fit reveals that the decay 
constant 72 of the second mode is about a factor of two larger as compared to the first mode while the obtained value for u 
differs rather slightly from the first mode's one. The light blue line is left after subtracting both modes from the FID and 
reveals that within experimental resolution the FID is accurately described by /i + /2 for times greater than 75 lis or for a 
decay over four orders of magnitude. 



150 jis, C.f. grey curve in Fig. [lib). The red line in 
Fig. [TJb) is a fit according to Eq. (|2| based on our decay 
between 160 and 250 jis using the nonlinear Levenberg- 
Marquardt algorithm [121 ISO]- This window was cho- 
sen since we may expect a faster second mode and since 
the algorithm should not operate in low signal-to-noise 
ratio. We verified that a slight change in the starting 
point did not change the obtained fit parameters signif- 
icantly. The obtained values are given in Tab. |T] as FID 
(110) Ai. We find a significant deviation between our 
value (71 — 40.3 ms~^) and the ones determined from 
Engelsberg and Lowe's data, but also those by Sorte et 
al. (~ 43ms~^). We also note a substantial difference 
between the fit (red line) and the data for shorter times. 
This is shown more clearly in Fig. [2]ja) where the green 
curve is the difference between the actual decay and the 
fit to the long-time decay (first mode Ai = 71 -f iwi). 
Clearly, the green curve in Fig. [2ja) is strong evidence 
for Fine's second mode, since it is approximately expo- 
nential with a decay rate of roughly 271 (It disappears 
in the noise at about 140 jis while the first mode has dis- 
appeared at 300 jis). Note that this result could not be 
obtained by earlier measurements due to the lower signal- 
to-noise ratio that did not permit a fit at sufficiently long 
times. The significantly higher values for 7 obtained from 
earlier measurements are due to the influence of the sec- 
ond mode. 

The parameters of the second mode (green curve) in 
Fig. [2][a) are estimated from a fit to Eq. ( 2]) between 75 
and 135 ys, shown as red curve in Fig. ^^b). The blue 



curve in Fig. shows the residuum, revealing that the 
FID can be described accurately in terms of two modes 
for times greater than 70 jis. The estimates for the second 
mode are given as FID (110) A2 in Tab.|lj 

While we are convinced that the green curve in 
Fig. [SJa) is evidence for Fine's second mode, we decided 
on an independent check, i.e., the measurement of the 
FID of CaF2 for a different orientation of the crystal 
with respect to the applied magnetic field. The resulting 
data and analysis are summarized in Fig. |3] Due to the 
dependendence of the secular part of the dipolar Hamil- 
tonian on the orientation of the applied magnetic field, 
a different shapes of the FID are obtained for different 
orientations. For the (100) orientation the FID decays 
faster and oscillates at a higher frequency. The data can 
again be understood in terms of two modes with the de- 
cay rate of the second mode being about two times larger 



Measurement 7 (ms 



(rad ms 



Engelsberg 


43.0(1) 


101.9(5) 


Sorte 


43.3(1) 


106.2(1) 


FID (110) Al 


40.3(2) 


103.5(2) 


FID (110) A2 


92(3) 


85(1) 


FID (100) Al 


48.4(3) 


154.0(2) 


FID (100) A2 


109(2) 


142(2) 



TABLE I. Obtained values for 7 and oj. All values were ob- 
tained using a non-linear least squares approximation. For the 
FID at 20 K deviations from the fit are due to thermal noise 
(c.f. Fig[2|, hence 95% confidence intervals were estimated 
from the covariance matrix. 
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FIG. 3. The same plot as in Fig. [2] but for B parallel to (100). The FID (a) is again measured with different attenuations; the 
first 23 lis are measured using a spin locking sequence, the time intervalls 24 to 51 11s and 52 to lllps are measured with 25 
and 12 dB, respectively. A Gaussian line narrowing of 426 jis was applied to remove the influence of the field inhomogeneity. 
The fit of the first mode /i is performed in the time intervall 100 to 200 jis. The residue, shown as green line in (a), reveals 
that in this time interval the data is well described by a single mode and again shows oscillatory decaying behavior for earlier 
times. This oscillatory behavior can well be described in terms of a second mode /2, plotted in (b), with a decay rate about 
two times the first mode's one. 



than the rate of the first mode. The fit results for both 
modes are given in Tab. |lj 

To conclude, we have studied the long-time behavior of 
a macroscopic system of dipole coupled spins 1/2. The 
data is interpreted in terms of Fine's theory based on 
the notion of microscopic chaos. While earlier work is 
in agreement with this theory, so far only a single decay 
mode was found. Therefore, it is difficult to discriminate 
other theories using a memory function approach. By 
increasing the signal-to-noise ratio by two orders of mag- 
nitude, compared to so far available measurements, we 
can determine the first decay mode with high accuracy. 
After subtracting this mode from the FID we are left 
with a second universal decay mode that decays about 
two times faster than the first mode. Our findings thus 
support Fine's theory that predicted a well-isolated sec- 
ond decay mode and correctly estimated the difference 
between the decay rates of the first and the second mode. 
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